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ABSTRACT 

We present results for the estimation of gravity, effective temperature and radial velocity of poorly 
studied chemically peculiar stars recently observed with the spectropolarimeter ESPaDOnS at CFHT 
in the frame of the VeSElkA (Vertical Stratihcation of Element Abundances) project. Eor four of the 
stars from our sample (HD 23878, HD 83373, HD 95608 and HD 164584), the values of their effective 
temperature and surface gravity are determined for the very first time. Grids of stellar atmosphere models 
with the corresponding fluxes have been calculated using version 15 of the PHOENIX code for effective 
temperatures in the range of 5000K to 15000K, for the logarithm of surface gravities in the range of 
3.0 to 4.5 and for the metallicities from -1.0 to h- 1.5. We have used these fluxes to fit the Balmer line 
profiles employing the code EITSB2 that produces estimates of the effective temperature, gravity and 
radial velocity for each star. When possible, our results are compared to those previously published. 
The physical characteristics of 16 program stars are discussed with the future aim to study the abundance 
anomalies of chemical species and the possible vertical abundance stratification in their stellar atmosphere. 

Subject headings: atomic processes - line: profiles - stars: abundances - stars: atmospheres - stars: chemically 
peculiar 


1. Introduction 


Despite of the fact that some chemically peculiar 
(CP) stars have ’’stable” atmospheres, they sometimes 
also show variability of their spectra with the period 
of stellar rotation due to the horizontal inhomoge¬ 
neous distributi on of elements abu ndance in their stel¬ 
lar atmosphere ( Khokhlo'^l 1975 ). A significant frac¬ 
tion of CP stars shows sign atures of strong magnetic 
fields ( Bychkov et al.l 2003 ) and their structure cor¬ 
relates with the patches of ov erabundance (or under¬ 
abundance) of some elements ( Kochukhov et al.ll2002 ; 
Shavrina et al .I2OIOI) . It appears that the magnetic field 
can intensify accumulation or de pletion of chemical 


elements at certain optical d epths (IRvabchikova et al 
2008; Alecian & Stifl 2010 ). Some stars also show 


'Based on observations obtained at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research 
Council of Canada, the Institut National des Sciences de TUnivers 
of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii. 


vertical stratification of abundance of several chemical 
speci es ( Savanov & Kochukhov! 1998 ; Rvabchikova et al 


2004t iKhalack et al.ll2013h . which can be ex plained in 
terms of the mechanism of atomic diffusion ([Michaud 
19701) . 


Accumulation or depletion of chemical elements at 
certain optical depths brought about by atomic diffu¬ 
sion can modify the structure of stellar atmospheres 
(Hui-Bon-Hoa et al. 2000; LeBlanc et al. 2009, 2010) 
and it is, therefore, important to gauge the intensity 
of such stratification. The slowly rotating CP stars 
with V sini < 40 km s^^ are good candidates to study 
the abundance stratification of elements with optical 
depth. Their small rotational velocities result in com¬ 
paratively narrow and unblended line profiles which 
are well suitable for abundance analysis. Based on 
their slow rotation, we may assume the hypotheses of 
a hydrodynamic ally stable atmosphere in these stars 
which is necessary for the diffusion process to take 
place. 

The diffusion velocities of different chemical 
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Table 1: List of the observed slowly rotating CP stars. 


Star 

m\ 

At 

(s) 

S/N 

HD 15385 

6.2 

1600 

1100 

HD 22920 

5.5 

920 

1150 

HD 23878 

5.2 

660 

950 

HD 53929 

6.1 

1160 

1000 

HD 68351 

5.6 

920 

1100 

HD 71030 

6.1 

1140 

1100 

HD 83373 

6.4 

1524 

1000 

HD 90277 

4.7 

520 

1000 

HD 95608 

4.4 

416 

1300 

HD 97633 

3.3 

152 

1300 

HD 110380 

3.6 

200 

1300 

HD 116235 

5.9 

1040 

870 

HD 164584 

5.4 

880 

1300 

HD 186568 

6.0 

1300 

1000 

HD 209459 

5.8 

1160 

1000 

HD 223640 

5.2 

680 

1200 


species depend on the relative values of gravity and 
radiative acceleration resulting from the momentum 
transfer from the radiation field to these chemical 
species (e.g. Gonzalez et al. 1995). The momen¬ 
tum transfer depends on the opacity of the species 
under consideration and on the local monochromatic 
radiation field, which, in turn, via the monochromatic 
opacities, depends on the local abundances of the dif¬ 
ferent species. In a hydrodynamically stable atmo¬ 
sphere, the diffusion process may result in vertical 
stratification of chemical abundances. Detection of 
vertical abundance stratification of chemical species 
in atmospheres of slowly rotating CP stars is an indi¬ 
cator of the effectiveness of the diffusion mechanism 
responsible for the observed peculiarities of chemical 
abundances. Comparing the observed characteristics 
of vertical stratification of chemical abundances with 
the results of theoretical modeling will help to verify 
and improve the self-consistent models of stellar atmo¬ 
spheres. Such models were calc ulated with a modifie d 
version of the PHOENIX code dLeBIanc et al. 2009h . 
These models were mostly applied to blue horizontal- 
branch stars ( LeBlanc et al.l 2010l) . Stift & Alecian 
(120 12h also have calculated chemically stratified atmo¬ 
spheric models using the CARATSTRAT code applied 
to ApBp stars. 

Vertical stratification of the chemical abundances in 
stars can be estimated through the analysis of multiple 


lines that belong to the same ion of the studied ele- 
ment dKhalack & Wadell2006HKhalack et al.ll2007l) i 


ing the ZEEMAN2 code ( Landstreed 1988 ). This pro¬ 
cedure was successfully implemented to detect verti¬ 
cal abundance stratification in the atmospheres of sev¬ 
eral blue horizontal-branch stars (Khalack et al. 2007, 
2008, 2010). 


Therefore, we have initiated a new project enti¬ 
tled ’’Vertical Stratification of Elements Abundances” 
(VeSElkA - meaning rainbow in Ukrainian) aimed to 
search for, and study the signatures of abundance strat¬ 
ification of chemical species with optical depth in the 
atmospheres of slowly rotating CP stars. A list of 
suitable candidates has been compiled based on the 
catalog of Ap, HgMn and Am stars of ?. Prom the 
very beginning, we have concentrated our attention 
on the relatively bright, slowly rotating and poorly 
studied CP stars that can be easily observed with ES- 
PaDOnS (CPHT) in spectropolarimetric mode. Up 
to now, we have obtained high signal-to-noise ratio 
and high resolution spectra for a sample of selected 
CP stars (see Sec. EJ and found indication of verti¬ 
cal stratificat ion of iron abundan ce in HD 95608 and 
HD 116235 ( Khalack et al. 2013 ). Not all stars in our 
sample are expected to contain vertical stratification 
of elements in their atmosphere. Por instance, AmPm 
stars are expected to have chemically h omogeneous 


atmospheres due to conv ective mixing (iRicher et al 


2000t iRichard et al.ll2001h . Our analysis will however 


give estimates for the average abundances in the atmo¬ 
spheres of these stars. Such results could help confirm 
or disprove their classification as CP stars. 


This paper aims to present results for the determi¬ 
nation of the fundamental parameters (Teff and logg) 
as well as radial velocity of 16 selected stars. This is 
the first estimation of Tes and logg for four of them. 
The observations and the reduction procedure are de¬ 
scribed in Section |3 The models used and the fitting 
procedure are discussed in Sections [3] and |4] respec¬ 
tively. The main results are presented in Section |3 to¬ 
gether with the description of the properties of each 
studied star. A discussion follows in Section|6l 


2. Observations 

High resolution (R=65000) Stokes IV spectra of 
several CP stars with Vsin/ < 40 km s^* were ob¬ 
tained recently with ESPaDOnS (Echelle SpectroPo- 
larimetric Device for Observations of Stars) employ¬ 
ing the deep-depletion e2v device Olapa. ESPaDOnS 
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Fig. 1.— The observed spectra (thick line) of 108 Vir (left) and a Leo (right) are well fitted by synthetic spectra (thin 
dotted line) that corresponds to T^s = 9715 K, logg = 3.56, [M/H]= -0.5 (X^/v = 2.64) and = 11890 K, logg = 
3.68, [M/H]= 0.0 (xVv = 12.26) respectively. 


allows the acquisition of an essentially continuous 
spectrum throughout the spectral range from 3700A 
to 10500A in a single exposure ( Donati et al. 19991) . 
The optical characteristics of the spectrograph as 
we ll as the ins t rumen t performances are described 
by Donati et al. ( 2006[ \ The obtained spectra were 
reduced using the ded i cated software package Libre- 
ESpRIT ( Donati et al. 1997 ) which yields both the 
Stokes I spectrum and the Stokes V circular polari¬ 
sation spectrum. To infer the effective temperature 
and gravity of the observed stars, we have used their 
non-normalized spectra. 

Table [Upresents a list of the observed slowly rotat¬ 
ing chemically peculiar stars. The hrst and the second 


^ For more details about this instrument, the reader is invited to visit 
www.cfht.hawaii.edu/Instruments/Spectroscopy/Espadons/ 


columns provide respectively the name of the star and 
its apparent visual magnitude, while the third and the 
forth columns contain the accumulation time and the 
maximal signal-to-noise ratio for the acquired Stokes I 
spectra. 


3. Grid of models 


A new library of high resolution synthetic spec¬ 
tra has been created in order to determine the effec¬ 
tive temperature and gravity of stars observed in the 
frame of the project VeSElkA. A grid of stellar at¬ 
mosphere models and corresponding fluxes have been 
calculated using versio n 15 of the PHOENIX code 
( Hauschildt et al. 1997h for 5000 K ^ Teff ^ 15000 
K and 3.0 ^ logg ^ 4.5. Eor the effective tempera¬ 
tures from 5000 K to 9000 K, we have used a 250 K 
step, while for the higher temperatures up to 15000 K 
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a 500 K step was used. The gravities have been cal¬ 
culated with a step 0.5. The theoretical fluxes have a 
resolution of 0.05A in the visible range from 3700A 
to 7700A. We have calculated grids of models for so¬ 
lar metallicity ( Grevesse et al.ll201^ as well as for the 
metallicities [M/H]= -1.0, -0.5, h-0.5, h-1.0, h-1.5. For 
all the grids, the microturbulent velocity is assumed to 
be zero. The synthetic spectra have been corr ected for 
the air wavelength using the equation given by MortonI 
(Il99lh . 

Exploiting the possib ilities of the PHOENIX code 


(iHauschildt et al.l Il997h . we have calculated spheri¬ 


cally symmetric models of stellar atmospheres for 
which the gravity depends on the stellar radius. In 
this case, all the structure of the atmosphere can be 
calculated knowing the effective temperature the 
surface gravity logg and the stellar mass M*. Eor the 
given values of effective temperature and surface grav¬ 
ity, we have derived the respective stellar mass trough 
the interpolation of data fo r physical pa rameters of 
normal main-sequence stars (|PoppeJl98d) . 


To verify the accuracy of the grids of stellar at- 


ESPaDOnS in the spectropolarimetric mode. The pro¬ 
files of nine Balmer lines in the non-normalised spec¬ 
tra of a Leo and 108 Vir have been fitted with the 


help of the EITSB2 code (iNapiwotzki et al.ll2004f) us¬ 
ing the grid of models calculated for the metallicity 
[M/H]=-0.5 and 0.0 respectively (see Eig.[T]). In order 
to model a continuum in the vicinity of the analysed 
Balmer lines, the code EITSB2 uses a linear function 
fitted over the line profile region and the continuum as 
a part of the fitting procedure. This approach provides 
a better fit quality than the standard procedure in which 
one tries to fit the continuum regions (relatively free of 
lines) to the left and to the right of the anal ysed Balmer 
line profile separately ( Naniwotzki Il20 14 ). 

I n the case of a L eo, assuming its V sin / = 300 km 


(lAto et alJl20()2 l). the best fit is obtained for solar 


metallicity, Teff = 11890 K and logg = 3.68, result¬ 
ing in X^/v = 12.26 (see right panel of the Eig.[T]). A 
similar simulation that employs th e grid of models cal- 
culated with PHOENIX-16 code ( Husser et al.ll2013li 
gives the best fit for the solar metallicity. Teg = 11976 
K and logg = 3.75 (X^/v = 11.57). These results are 
close to each other and to the values T^g = 11962 ± 


^Spectra of 108 Vir and of a Leo where downloaded from the CADC 
database via http://www4.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/en/cfht/ 


115 K and logg = 3.56 obtained by Gray et aP (l2003h 
from fitting low resolution spectra of a Leo and its 
photometric data. Meanw hile, using the Atlas9 grid^ 
( Castelli & Kuruczl2004 ) to fit the Balmer line profiles 
in the same spectrum of a Leo, the best fit is obtained 
for the metallicity [M/H] = -0.3, Teff = 12660 K and 
logg =3.80 (X^/v = 13.76). It appears that for a Leo 
the use of the Atlas9 grids leads to the best fit with 
comparatively higher value of and to higher val¬ 
ues of effective temperature and surface gravity. 

In the case of 108 Vir, assuming its Vsin/ = 83 


km s ' dAmmler-von Eiff & Reiners 2012 ). the best 


fit is obtained for the metallicity [M/H]=-0.5, Teff = 
9715 K and logg = 3.56, resulting in X^/v = 2.64 
(see left panel of Lig. [T]i. The grids calculated with 
PHOENIX-16 code and the Atlas9 grids result in the 
best fits with [M/H] = -0.2, T^g = 10280 K, logg = 
3.86 (xVv = 2.57) and [M/H] = -0.3, T^g = 9760 
K, logg = 3.64 (x^/v = 2.88) respectively. Our es¬ 
timate of the effective temperature is close the one 
obtained wit h the Atlas9 grids and to 7 ].ff = 9840 


K found by lAmmler-von Eiff & ReinersI (12012h . but 


spectra of a Leo 

is significantly smaller than the effective temperature 

jj obtained with 

obtained using the PHOENIX-16 grids ( 

Husser et al. 


2013h . Our estimate of the surface gravity is similar to 


the one obtained for the Atlas9 grids but smaller than 
the one obtained with the PHOENIX-16 grids. 

By comparing our results for a Leo and 108 Vir 
with the relatively well-established values for Tgg and 
logg, we may roughly estimate the uncertainties of the 
results presented here as ±200 K for effective temper¬ 
ature and ±0.2 for surface gravity. 

4. Fitting procedure 

4.1. Sensitivity to the set of analysed Balmer lines 

In order to study the sensitivity of the determined 
values of Tgg and logg to the set of analysed Balmer 
line profiles, we have compared the best fit results ob¬ 
tained when one of the Balmer lines is omitted from 
the analysis with the best fit results obtained from the 
analysis of all nine Balmer line profiles for HD 129956 
(108 Vir) and HD 87901 (a Leo) using the Phoenix- 
15 grids and the Atlas9 grids. The fitting proce¬ 
dure have been performed wi th the help of PITSB2 
code ( Naniwotzki et akl 20041) employing the afore¬ 
mentioned grids of stellar atmosphere models together 


^The Atlas9 grids are available at the 
http://zuserver2.star.ucl.ac.uk/~idh/NewGrids/Atlas9.C04/ 


4 


















































Fig. 2.— Variation of the best fit results obtained, when one of the Balmer lines (indicated at the X-axis) is omitted 
from the analysis, with respect to the values of Teff (left) and logg (right) obtained from the analysis of all 9 Balmer 
line profiles (see Fig.[T]l for HD 129956 (108 Vir) using the Phoenix-15 grids (open circles) and the Atlas9 grids (filled 
circles), and for HD 87901 (a Leo) using the Phoenix-15 grids (open diamonds) and the Atlas9 grids (filled diamonds). 
The error bars represent the internal errors of the fitting procedure with the F1TSB2 code and are much smaller than 
the true errors (see Section^. 


with the respective simulated spectra. The left panel 
of Fig. |2] shows the differences between the values of 
Teff obtained by omitting one of the Balmer lines from 
the analysis and the value of Teff given in Section |3] 
for HD 129956 and HD 87901 respectively. The cor¬ 
responding differences for logg are shown in the right 
panel of Fig.|2]for the same stars. In the case of Teff and 
logg, the zero difference corresponds to the respective 
values obtained using the Phoenix-15 grids and pre¬ 
sented in Fig. [U for each reference star. Meanwhile, 
the results shown in the Fig.|2]for the Atlas9 grids are 
compared to the corresponding values of stellar param¬ 
eters found for the reference stars in Section |3] using 
the Atlas9 grids. 

For HD 129956 (108 Vir), we can see that its best fit 
values of Teff and logg jump up if one uses the Atlas9 
grids and excludes from the analysis one of the follow¬ 
ing Balmer lines H^, H^, 7/8 or H9. Meanwhile, the 
use of the Phoenix-15 grids results in relatively small 
variation of Teff and logg when one of these Balmer 
lines is excluded from the analysis. 

In the case of HD 87901 (a Leo), application of 
the Atlas9 grids or of the Phoenix-15 grids leads, un¬ 
der the same conditions, to relatively small variations 
of Teff and logg that do not exceed 100 K and 0.06 
respectively (see Fig. EJ. The errors chosen for Teff 
(±200 K) and logg (±0.2) therefore seem reasonable. 
It appears that the use of the Phoenix-15 grids results in 


more stable best fit data for Teff and logg with respect 
to the set of Balmer line profiles employed for the anal¬ 
ysis. Meanwhile, the use of the Atlas9 grids may cause 
significant errors in the determination of fundamental 
stellar parameters if the effective temperature is close 
to 10000 K, depending on which set of Balmer lines is 
used. 

The variations of Teff and logg are most sensitive to 
the Hp, Hy, 7 / 5 , and T/g Balmer lines when one employs 
the Phoenix-15 grids. It seems that this sensitivity does 
not vary much in the range of effective temperatures 
from 9700 K to 12000 K (see Fig.Ejl. 

4.2. The best fit results 

Nine Balmer line profiles have been fitted in the ob¬ 
served spectra of the selected slowly rotating CP stars 
to find their effective temperature and surface gravity. 
For each star, the fitting procedure has been performed 
for the metallicities [M/H]= -1.0, -0.5, 0.0, h-0.5, h-LO 
using the Phoenix-15 grids. Among the obtained re¬ 
sults, the fundamental parameters associated to the fit 
with the smallest value are chosen. These best 
fit results for effective temperature, surface gravity, ra¬ 
dial velocity and metallicity are presented respectively 
in the second, third, fourth and fifth columns of Ta¬ 
ble |2] together with the fit quality in the sixth column. 
The values of the effective temperature and the surface 
gravity found for these stars by other authors are given 
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Table 2: Determined values of the effective temperature and the surface gravity for the programm CP stars. 


Star 

Tss 

(K) 

Balmer lines 
logg Vr 

(kms-') 

[M/H] 

xVv 

Preyious publications 

Teff logg Vsini 

(K) (kms-i) 

HD 15385 

8230±200 

4.00±0.2 

22.0±1.0 

0.0 

0.65 

8154“ 

4.12“ 

21“, 29" 

HD 22920 

13640±200 

3.74±0.2 

20.0±2.0 

-0.5 

5.67 

13700“ 

3.72“ 

30“', 39" 

HD 23878 

8740±200 

3.86±0.2 

29.5±1.0 

0.0 

0.62 



24" 

HD 53929 

13950±200 

3.90±0.2 

15.0±1.0 

-1.0 

3.20 

14050±250“ 

3.60±0.25“ 

25", 30“ 

HD 68351 

10080±200 

3.22±0.2 

18.1±1.0 

0.0 

1.17 

10290±340^ 


33" 

HD 71030 

6780±200 

4.04±0.2 

38.1±1.0^* 

0.0 

0.28 

6541±47S 

4.03±0.05® 

9±2" 

HD 83373 

9800±200 

3.81±0.2 

26.5±1.0 

0.0 

0.92 



28" 

HD 90277 

7250±200 

3.62±0.2 

14.5±1.0 

0.0 

1.29 

7440' 

3.46' 

34" 

HD 95608 

9200±200 

4.25±0.2 

-10.4±1.0'' 

h-0.5 

0.59 



2l",17±2" 

HD 97633 

8750±200 

3.45±0.2 

8.2±1.0 

0.0 

0.61 

8790±35F' 

3.59±0.89^' 

23" 

HD 110380 

6980±200 

4.19±0.2 

-17.6±1.0 

0.0 

0.31 

6720^ 

4.20^^ 

23" 

HD 116235 

8900±200 

4.33±0.2 

-10.3±1.0'' 

h-0.5 

0.49 

8570' 

4.23' 

20±2" 

HD 164584 

6800±200 

3.54±0.2 

-11.2±1.0 

0.0 

1.16 




HD 186568 

11070±200 

3.44±0.2 

-9.5±1.0 

-0.5 

1.81 

11596±120” 

3.39±0.15” 

18” 

HD 209459 

10310±200 

3.62±0.2 

-0.3±1.0 

0.0 

0.97 

10455 ±400” 

3.52±0.15” 

14" 

HD 223640 

12500±200 

4.08±0.2 

17.0±2.0 

H-1.0 

1.65 

12429±435« 

3.93±0.23« 

28" 


Notes: fkiinzli & Northl llQ^. jRover et alJI2002^■ ICatanzaro et alJll999l). iLeone & Manfr^tl996ft. ISmlth & Dworetskil 

IPrugniel et aljf20Tlh . iKhalack et alh2013l) . fBerthell IT99QI) . j^leva_&_N^izdeM^j^l2l). IBoesgaard & Trvpiccolll98tll . Erspamer & NorthI t20Q3l) . 

" iHubrig & CastelliuioOll) 


in the seventh and eighth columns respectively. Mean¬ 
while, the previously published values for V sin / are 
presented in the ninth column. Examples of best fits 
for HD 23878 and HD 164584, and for HD 95608 and 
HD 209459 are shown in Fig.|3]and Fig. Irrespectively. 


5. Individual stars 


5.1. HD15385 


Kunzli & NorthI (Il998h have classified HD 15385 
(HR 723) as a A2m star having a mass M* = 1.85M0 
and age logf = 8.79±0.12. Taking into account its rela- 
tivel y small value of rot ational velocity V sin i = 29 km 
s^* dRover et al.l l2002l) . this star was chosen for the 
VeSElkA project. Our estimate of the effective tem¬ 
perature Teff = 8230±200 K and the surface gravity 
logg = 4. 0ifc0.2 are in good agre ement with the values 
found by Kunzli & North ( 19981) (see Tab.|2]l- Our es¬ 
timate of the radial velocity for HD 15385 (see Tab.|2ll 
is similar to the value of Vr = 21.3 km s^' obtained 


by IWilsonI (Il953h . but is signific antly higher than th e 
value Vr = 15 km s^* reported bv IPalmer et al.l (Il968l) . 


5.2. HD22920 


The silicon star HD 22920 (22 Eri, HR 1121) shows 
a wea k photometric va riability with the period P = 
3‘^.95 ( Bartholdy 1988 ) and the presence of a rather 
weak longitudinal m ag netic field Ri = 3J_0±160 G 
( Bychkoy et al. 2003 ). Catanzaro et al. (1999) haye 
found Teff = 13700 K and logg = 3.72 for this star. 
Similar parameters (rp jf = 13700 K and logg = 3.72) 
haye been adopted by iFeone & Manfrd ( 19961) to re¬ 
produce the obseryed spectrum of HD 22920. Our esti¬ 
mate of the effectiye temperature and the surface gray- 
ity (see Tab.|2]i are consistent with these results. Mean¬ 
while, the radial yelocity Vj- = 20.0±2.0 km s^* ob¬ 
tained in the present study seems to be higher than the 
preyiously reported yalues 16.4±3.1 km s^' (Feone 
& Manfr e 1996; Widen et al. 2000) and 17.2±0.7 
km s^' ( Gontcharoy 200^. From the a nalysis of 
HD 22920 spectra, Feone & Manfre ( 19961) haye also 
found its rotational yelocity V sini = 30 km s^^ 


We haye acquired seyeral Stokes IV spectra for this 
star that show a strong yariability of Si II, Ti II, Cr II, 
Fe II line profiles with rotational phase, while Mg II 
and He I line profiles seem to be less yariable. Our high 
S/N spectra (see Tab. [T]i reyeal a weak yariability of 
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Fig. 3.— The observed spectrum (thick line) of HD 164584 (left) and HD 23878 (right) is well fitted by synthetic 
spectrum (thin dotted line) that corresponds to Tes = 6800 K, logg = 3.54 (x^/v = 1.16) and Teg = 8740 K, logg = 3.86 
(X^/v = 0.62) respectively. 



the He I 5876A line profile for which ICatanzaro et al 
(Il999h have found no variability. 


value Vr = 28.4±0.5 km s * obtained by iGontcharov 
(l2006ll (see Tab.l2Tl. 


5.3. HD 23878 

Based on the measurement of the anoma lously 
low line strength ratio Sc II 4246A/Sc II 4215 A, Conti 
( 1965 ) has suggested that HD 23878 might be an Am 
star. HD 23878 (28 Eri, HR 1181) appears to be a 
low a mplitude variable s tar with a period of about 
7‘^.17 ( Mathvs et al.lll98^ having a rota tional velocity 
Vsini = 24 km s^' ( Rover et al. 2002). Meanwhile, 
Abt & Morrell ( 1995 ) have reported V sin / = 18 km 


s^^ for this star. HD 23878 has not previously been 
studied spectroscopically in detail and the fitting of 
the Balmer line profiles (see right panel at the Fig. [3) 
results in T^s = 8740±200 K and logg = 3.86±0.20. 
Our estimate of radial velocity is in agreement with the 


5.4. HD 53929 

HD 53929 is a HgMn star with mild but significant 
enhancement o f manganese abundance lo gAivin/A^H 
= -5.85±0.20 ( Smith & Dworetskv 1993h and with 
a mercury abundance \os Nut>/Nu = -9 .90±0.20 at 


the level of normal stars (iSmithI Il997l) . This star 


has also much lower abundance of chromium, cobalt 
and nickel as compared to their solar abundances 
( Smith & Dworetskv 19931) . The analysis of Balmer 
line profiles in the single spectrum available for 
HD 53929 results in an effective temperature which 
is in good agreement with t he temperature obtained by 
Smith & Dworetskv ( 1993 ) (see Tab. ID). Meanwhile, 
our estimate of the surface gravity logg = 3.90±0.20 
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seems to be a little higher than the value log ;? = 
3.60±0.25 reportedbv Smith & Dworetskv ( 1993 ). 
Rover et alJ (1200211 have obtained V sin i = 25 km 


5.6. HD 71030 


s ^ for HD 53929. Its radial velocity increased over 
the last 50 years going from 6.1 k m s~* (Evans 1967 ; 
Hube 1970), to 11.2±1.3 km s^^ (lGontcharovll2006h 


and to 15.Oil.0 km s^^ found in this study. Is seems 
that HD 53929 may be a member of a long periodic 
double system. 

5.5. HD 68351 


In the catalog of Renson & Manfroi^ (2009), this 
object is classified as a star of spectral class AO 
wi th the strong l ines o f Si and Cr. According 


to 


Abt & Morrelll (119951) . the rotational velocity of 
HD 68351 is Vsin; = 25 km s~*, while Rover et al 


(120021) and lAuriere et akl (l2007h have reported V sin / 
= 33 km s^*. Its rotati o nal pe riod P = 4"^. 116 was 
determined by IStepieh (119681) . In this study, we 
present the results for its fundamental parameters 
= 10080i200 K and logg = 3.22i0.20 based on the 
analysis of Balmer line profiles which ar e in agreement 


with t he data previously published by lAuriere et al. 


(l2007l) for this star. Our estimate of the radial velocity 
seems to be a lit tle smaller tha n the value 14 = 19.9 km 
s^^ obtained by Evans ( 1967 ) (see Tab.|2]). 

Taking into account the high values of luminos¬ 
ity L* = 466 ifc 225LfTi and radius —6.6 iL2.4Rq 
of HD 68351 ( Auriere et al. 20071) and its low sur¬ 
face gravity logg = 3.22±0.20, this star is beyond the 
end of the main sequence a nd most prob a bly be longs 
to the III luminosity class. lAuriere et al.l (l2007l) have 
detected the presence of a longitudinal magnetic field 
B/ = 325 ±45 G from the analysis of HD 68351 LSD 
profiles. 

For this star, we have obtained two Stokes IV spec¬ 
tra with a time gap of 2 days. The Stokes I spectrum 
of HD 68351 appears to be strongly variable, while 
the Stokes V does not show signatures of the pres¬ 
ence of a stro ng magnetic field. Th is is consistent with 
the results of lAuriere et al.l (l2007h that have observed 
HD 68351 15 times, but have found the magnetic field 
signatures in only 5 spectra. A preliminary analysis 
of the two obtained spectra results in slightly differ¬ 
ent values of radial velocity 14= 18.1±1.0 km s^ ^ and 
14.7±1.0 km s^* that can be an indicator of possible 
binarity of HD 68351. 


Cenarro et al.l (l2007l) have reported HD 71030 
(25 Cnc, HR 3299) to be a main-sequence star of spec¬ 
tral class F6. Based on th e analy sis of a low resolution 
spectrum, [BalachandranI ( 1990h has reported higher 
than solar lithium abundance logNu/Nn = -9.37 and 
slightly underabundant iro n logVpe /A4t = -4.80 ±0.05 
for this star. Meanwhile, Khalack et al.l (l2013h have 
found nearly solar abundance of Fe II logNpe/Na = 
-4.43±0.34, Cm logNa/Nn = -6.27±0.09 and Nil 
logVNi/A4i = -5.77±0.26. They have also found that 
Fe, Cr and Ni are uniformly distributed with respect to 
optical depth in the atmosphere of this star. 

Our estimation of the surface gravity for HD 71030 
is the same as the value logg = 4.03±0.05 published 
by IPrugniel et alJ (1201 ll) (see Tab. IH, while our esti¬ 
mation of the effective temperature Teff = 6780±200 K 
se ems to be sligh t ly high er as compared to the results 
of IPrugniel et al.l (|201 ih (but within the given error 
bars). The value for the radial velocity determined 
in this work is in agreement with the previously ob- 
tained results of 37.4 k m s~^ dN ordstro rn et al. 2004) 
and 37.1±0.4 km s^^ ( Gontcha rovI l2006 l) . T he esti¬ 
mate of V sin / = 9±2 km s~* dKhalack etak 2013 ) is 
also in good agreement with the value V sin / = 8 km 
s^^ provided by Balachandran (1990). 

5.7. HD 83373 


HD 83373 (34 Hya, H R 3832) is reported by 
Renson & Manfroid ( 20091) as a chemically peculiar 
star of spectral class Al that has an enhanced silicon 
abundance. Royer et al. (2002, 2007) have measured 
its rotational velocity Vsin/ = 28 km s^^ This star 
has not been previously studied spectroscopically and 
we provide here first estimates of its effective tem¬ 
perature Teff = 9800±200 K and surface gravity logg 
= 3.81 ±0.20. Our estimate of the radial velocity for 
HD 83373 (see Tab. |2]i is in good agree ment with the 
value 14 = 26.9±0.5 km s^' reported bv iGontcharovI 

(l2006l) . 


5.8. HD 90277 

Abundance anal ysis of HD 902 77 (HR 4090) has 
been performed by Berthe) ( 1990l) who has found a 
strong overabundance of Y, Zr and Ba, while Ti, Cr, 
Mn, Fe and Ni appear to be slightly overabundant. 
Royer et al. (2002, 2007) have reported V sin / = 34 
km s^* for HD 90211. Our estimates of the effective 
temperature and the surface gravity are similar to the 
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Fig. 4.— The observed spectrum (thick line) of HD 95608 (left) and HD 209459 (right) is well fitted by synthetic 
spectrum (thin dotted line) that corresponds to Teff = 9200 K, logg = 4.26 (X^/v = 0.59) and Teff = 10310 K, logg = 
3.62 (x^/v = 0.97) respectively. 


values obtained by Berthed (Il990ll taking into account 
the error bars (see Tab.|2]l. The value for the radial ve¬ 
locity determined in this work agrees well with the pre- 
viously obtained resu lt 14 = 13.7 km s~* ( Evans! 1967 ) 


and 13.7±0.6 km s ^ Gontcharovl ( 200^ . 


5.9. HD 95608 


Cowley et al.l (119651) have classified HD 95608 
(60 Leo, HR 4300) as a CP star of spectral class Aim. 
From the fitting of Balmer line profiles of HD 95608 
(see left panel of Fig. @1, we have found its effective 
temperature Teff = 9200±200 K and surface gravity 
logg = 4.25±0.20 (see Tab. |2|. The value for radial 
velocity determined in this work Vr = -10.4±1.0 km 
s^^ agrees we ll with the previou sly obtained results 
-10. 1 km s-^ dWielen et ai.ll2000l) and -11.1±0.7 km 
s^^ ( Gontcharovll2006 ). The rotational velocity V sin/ 


= 17.2±2.0 km s * obtained bv lKhalack et alJ (120131) 
for HD 95608 seems to be smaller than the value V sin / 


= 21 km s ' reported by Rover et al.l (12002 ), but higher 
than y sin / = 13 km s^* published bv lAbt & Morrell 


(Il995h . 

Preliminary abundance analysis has shown that Ti II 
is slightly underabundant, while Fe I, Fe II are over- 
abund ant in this star as com pared to their solar abun¬ 
dance ( Khalack et al.ll201~3 ). This same study has also 
found that iron abundance appears to be vertically 
stratified in the stellar atmosphere of HD 95608. 

5.10. HD 97633 


Renson & ManfroidI(120091) have classified HD 97633 


(0 Leo, HR 4359) as a CP star of spectr al class A 2 
that has enhanced Sr and Eu abundances (lHilllll995l) . 
Rover et al. (2002) have found V sin / = 23 km s^* for 
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this star. Our estimates of its effective temperature 
and surface gravity (see Tab. |2]) are in good agree¬ 
ment w ith the results obtained bv iKoleva & Vazdekis 
( 2012 ). The radial velocity Vr = 8.2±1.0km s^' found 
in this work is consistent with the value Vy = 13 km 


reported by Gontchar^ ( 2006 ). 


5.11. HD 110380 

HD 110380 (HR 4826) is in a binary system with 
HD 110379 with an orbital period P=171.4yr and 
a separation a=3.15". HD 110380 has Vsini = 23 
km s^* ( Rover et al. 2002 |) and an enhan ced lithium 
abundance ( Boesgaard & Trvpicco 19861) . Our esti¬ 
mate of its effective temperature seems to be a lit- 
tle higher than the value Tj-tr = 6720 K obtained by 
Boesgaard & Trvpiccol ( 198 ct . but the surface gravity 
is in good agreement with the value logg = 4.2 they 
found (see Tab. |2]i. The value for the radial velocity 
determined in this work Vy = -17.6±1.0km s^* is con- 
sistent with the value V ,- = -19.5 km s^* obtained by 
Nordstrom et al. (2004). 


HD 164584 (7 Sgr, HR 6724) has a spectral class A6- 
F4. From the analysis of Balmer line profiles (see left 
panel of Fig. [3), we have obtained Teff = 6800±200 K 
and logg = 3.54±0.20. 

Our estimate of radial velocity (see Tab. |2]) is in 
good agree ment with the value K- = -11.2±1.6 km s^* 
obtained bv iGontcharovI (l2006h for this star. 


5.14. HD 186568 

HD 186568 ( HR 7512) belongs to the normal B- 
type stars (B9 II) ( Hubrig & CastellilboOl ). Neverthe¬ 
less, it was included by ? to the list of chemically pe¬ 
culiar stars. Fro m the analysis of equivale nt widths of 
iron line profiles iHubrig & Castellil (1200 Ih have found 
its solar abundance in the atmosphere of this star. 

The effective temperature obtained for HD 186568 
in this study is smaller than the value derived by 


Hubrig & Castellil (120011) . while our value for the sur¬ 


face gravity is almost the same as the one determined 
by these authors (see Tab. |2|. Our estimate for radial 
velocity is c l ose to 14 = -8.8±0.1 km s^^ reported by 


5.12. HD 116235 


Gontcharov 


(l2006f) and to t he val ue -8.0±1.4 km s 


Morrell & Aba (Il992h who have used this 


Cowley et al.l (119651) have classified HD 116235 


reported by 
star as a radial velocity standard. 


5.15. HD 209459 

HD 209459 (21 Peg, HR 8404) is a normal B-type 


(HR 5040, 64 Vir) as a CP star of spectral class 
A2m. The effective temperature Teff = 8900±200 K 
obtained in this study for HD 116235 appears to be 

higher than the previous ly published values 8570 K - sts^(B9.5 V) which is often used as a comparison star 

(|Erst3amer&Nort h|2003|) and 8373 K ( |Ammler-von Eiff ^ Budaill2000l) because of its shaiT-lined 

20121). Meanwhile, the surface grav ity obtained here is spectrum with V sin / = 4 km s-^ dSmithlM 9921) and ab- 

simito to logg = 4.23 reported by | Erspamer & Nor^H chemical abundance peculiarities. Based on 

( 2 003 |). The estim ate of F sin i = 20±2.0 km s equivalent widths of iron line profiles in 

(lKhalacketalJl20lJ is in a good agreement with the - . ^ spectrum of HD 209459. [Hubrig& Castellil (M 

yalue^l9.3±1.0kms_ ' obtained by| Ammler-von Eiff & Reinersf ^^^^ abundance that is close to its solar 

value. We also aim to use this object as a comparison 
star in the VeSElkA project. 

Our estimate of effective temperature and sur¬ 
face gravity (see Tab. |2|i obtained from the analysis 
of Balmer line profiles (see right panel of Eig. in 
is in agreement with the previously published re- 
sults Teff = 10455±400 K, logg = 3.52±0.15 of 


(12012h for this star. The radial velocity 14 = -10.3±1.0 
km s^^ obtained in this study is consiste d with the pre- 
viously published val ues -10.2 km s 
and -8.4±2.3 km s^^ ( Gontcharovl200i 



A preliminary abundance analysis of this star shows 
an en hanced abundance o f Ni I, Ee I, Ee II, Cr I and 
Cm (iKhalacke^^ 2013) in agr eement with the re¬ 
sults of [Erspamer & NorthI (2003), who have also re- 
porte d a strong overabund ance of Sr, Y and Ba. More¬ 


over, 


Khalack et alJ (1201 3h have also found signatures 


of vertical stratification of iron abundance in the atmo¬ 
sphere of HD 116235. 

5.13. HD 164584 


Hubrig & Castellil (1200 Ih . but smaller than the results 
T,ff = 11015±301 K , logg = 3.99±0.12 reported by 


Prugniel et al. (|201 ih . The radial velocity found in 


this study (see Tab.|2]) is also in a go od agreement with 
the value Vj- = 0 km s ^' reported by SmithI ( 1992 ). 


Renson & ManfroidI (120091) have reported that 
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5.16. HD 223640 


Renson & Manfroi3(l2009l) have reported HD 223640 


(108 Aqr, HR 9031) as a CP star of spectral class B9. It 
show s an overabundance of Si. Sr a nd Cr dNorth et al. 


1992h . iBailev & Landstreetl (12013h have found V sin/ 
= 31 ±3 km s^^ for this star and that Cr and Si abun¬ 
dances exceed their respective solar values more than 
by 1 dex. Analysis of polarization in Balmer lines re¬ 
veals the pres ence of a longitudina l magnetic field Bi 
= 643±218 G ( Bychkov et al J2003 ) which varies with 
period P = 3^ 73 dNorth et alJ 1 1992h . The same pe- 
riod P = 3 ^ 73523 9±0^.000024 has been found by 


North et alJ dl992h from photometric variability of 


HD 223640. 

The effective temperature and the surface gravity 
obtained in this study (see Tab. HI) are in good agree¬ 


ment with the results of iPrugniel et al.l (1201 Ih . but 


seem to be slightly different (but still inside the error 
bars) with respect to the values Tj-fp = 12300±500 K, 
logjg = 4.4±0.2 reported by IBailev & Landstre 
(12013). Our estimate of the radial velocity seems to 
be larger (but still within the error bars) than the valu e 
14= 12.7±2.8 km s^* reported by GontcharovI ( 200^ . 


6. Discussion 


Based on the catalogue of Ap, HgMn and Am stars 
(?) and the available measurement of rotational ve¬ 
locities of CP stars (Royer et al. 2002, 2007) we 
have compiled a list of CP stars (see Table [T]) suit¬ 
able for search of signatures of abundance stratifica¬ 
tion of chemical species with respect to optical depth 
in their atmospheres. These stars have been recently 
observed in the frame of VeSElkA project with the ES- 
PaDOnS at CFHT The method developed for the anal- 
ysis of vertical stratific ation of chemical abundance 
( Khalack & Wadell2006l) allows to determine the slope 
of abundance change relative to optical depth. 


In order to study the vertical stratification of chemi¬ 
cal abundance one needs to adopt an appropriate model 
for the stellar atmosphere and to determine its param¬ 
eters like effective temperature, surface gravity and 
metallicity. We were able to derive these parameters 
for the listed stars (see TablelU based on the analysis of 
nine Balmer lines (see, f or example. Fig. [3 and F ig. @1) 
using the FITSB2 code dNapiwotzki et ^ 20041) that 
employs a grid of theoretical fluxes calculated for dif¬ 
ferent values of Teff, logg and metallicity. Taking into 
account that the resolution R=65000 of the ESPaDOnS 
spectra, we have calculated a new library of synthetic 


spectra with a similar spectral resolution to properly 
fit the profiles of nine Balmer lines of stars observed 
in the frame of the project VeSElkA (see Section IH). 
Grids of stellar atmosphere models and corresponding 
fluxes have been calculated with the PHOENIX code 


dHauschildt et al.ll 19971) for 5000 K < T^ff 15000 K 


and 3.0 ^ logg ^ 4.5 for the metallicities [M/H]= -1.0, 
-0.5 ,h-0.5,h-1.0,h-1.5. 


We have used the spectra of 108 Vir (HD 129956) 
and a Leo (HD 87901) obtained with ESPaDOnS in 
the spectropolarimetric mode to verify the accuracy 
of our grids of synthetic fluxes. The results obtained 
for for these reference stars are close to the previously 
published data for their effective temperature and sur¬ 
face gravity (see Section^. 

In order to study the sensitivity of the determined 
values of Teg and logg to the set of Balmer lines used, 
we have compared the best fit results obtained when 
one of the Balmer lines is omitted from the analysis 
with the best fit results obtained from the analysis of all 
nine Balmer line profiles for the reference stars using 
the Phoenix-15 and Atlas9 grids. We have found that 
the use of the Atlas9 grids may produce some ambigu¬ 
ity in the determination of fundamental stellar param¬ 
eters if the effective temperature is close to 10000 K 
depending on which set of Balmer lines is used. Em¬ 
ploying our Phoenix-15 grids for simulations with dif¬ 
ferent sets of Balmer lines, we have shown that the 
estimates of Teff and logg are most sensitive to the //p, 
Ny, 7 / 5 , and T/g Balmer lines. This sensitivity does not 
change significantly in the range of effective tempera¬ 
tures from 9700 K to 12000 K (see Fig.|2]). 


We roughly estimate the uncertainties ±200 K and 
±0.2 respectively for the values of effective temper¬ 
ature and surface gravity. The relatively small vari¬ 
abilities found for and logg when using subsets 
of Balmer lines (see Subsection 14.11) . seem to show 
that this choice f or the errors i s reaso nable. It should 
be also noted that iHusser et al.l (l2013h have calculated 
grids of theoretical fluxes for Teg below 12000 K. 
Therefore we can not use those grids to determine Tgg 
and logg of all the stars selected for this study because 
some of them are much hotter than 12000 K (see Ta¬ 
ble |2]). Meanwhile, our grids of s ynthetic spectra pro- 
vide similar results as the grids of iHusser et al.l (1201 3l) 
for the selected stars with Tgg below 12000 K. Since 
our grids go up to Tgg = 15000 K, they can be applied 
for all of the stars in our sample. 


Our final results for effective temperature and sur¬ 
face gravity obtained for twelve of the stars presented 
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in Table|2]are consistent with the previously published 
data. For four other stars (HD 23878, HD 83373, 
HD 95608 and HD 164584), this study gives the first 
estimates of their and logg. 


It should be noted that some stars observed in the 
frame of the project VeSElkA have low Teff (see Ta- 
ble|2]) and won’t have vertical stratification of chemical 
species due to the presence of mixing due to convec¬ 
tion (Richer et al. 2000; Richard et al. 2001). Never¬ 
theless, they provide the possibility to verify the valid¬ 
ity of the method applied to analyse the vertical strat¬ 
ification of chemical abundances (IKhalack & Wade 


200611 and to determine their abund ance peculiarities. 
For instance, IKhalack et al.l (1201 3l) found no signa¬ 


tures of vertical stratification in HD 71030 and that 
the abundance of the analysed chemical species is 
close to their solar abundance. Our sample of anal¬ 
ysed stars includes the reference stars HD 186568 
and HD 209459 known to be a normal B-type stars 
dDworetsky & Buda) 2000ll with out abundance pecu¬ 


liarities (iHubrig & Castellill2001ll . These stars can also 
be used to verify the applied method and also to test the 
results for average abundances as well. The abundance 
analysis for all stars, even those where no stratification 
is found or expected, will be presented in upcoming 
papers. This analysis may be useful to confirm or 
disprove their CP type classification. 


Khalack et al.l (1201 3l) have already analysed several 


of the selected stars with the aim to detect vertical 
stratification of chemical abundances in their atmo¬ 
spheres and found signatures of vertical stratification 
of iron in HD 95608 and HD 116235. A detailed abun¬ 
dance analysis for the other stars of our sample is un¬ 
derway. We also plan to add more suitable stars to 
the project VeSElkA. A large amount of observational 
data relative to the vertical abundance stratification of 
chemical species obtained for CP stars with different 
effective temperatures will allows us to search for a de¬ 
pendence of the vertical abundance stratification rela¬ 
tive to the effective temperature, similar to the one that 
we have found for BHB stars (LeBlanc et al. 2010; 
KJialack et al. 2010). Our results will also be use¬ 
ful for comparison with theoretical modelling of verti¬ 
cal abundance stratification in stellar atmospheres (i.e. 
LeBlanc et al. 2009; Stift & Alecian 2012) and to bet¬ 
ter understand the diffusion process in the atmospheres 
of certain CP stars. 
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